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Abstract 

The nonperturbative parton distribution and wave functions are directly related 
to matrix elements of light-ray (nonlocal) operators. These operators are generaliza- 
tions of the standard local operators known from the operator product expansion. 
The renormalization group equation for these operators leads to evolution equations 
for more general distribution amplitudes which include the Altarelli-Parisi and the 
Brodsky-Lepage equations as special cases. It is possible to derive the Altarelli-Parisi 
kernel as a limiting case of the extended Brodsky-Lepage kernel. As new application of 
the operator product expansion the virtual Compton scattering near forward direction 
is considered. 



^ Supported in part by Deutschen Akademischen Austauschdienst and by Department of Energy contract 
DE-AC03-76SF00515. 



1 Light-Ray Operators in QCD 



It is generally accepted that nonlocal composite operators are necessary to describe liadrons 
in QCD. For example a meson operator M{xi,X2) could be represented by quark-antiquark 
fields connected by a path-ordered phase factor U{x2,Xi) : 

M(xi, X2) = HX2)TU (X2, Xi)^(Xi) (1) 

with 

U{x2,xi) = Pexp{-ig f ' dx^'A^ix)); (2) 

where F abbreviates 7- and A-matrices describing the spin and flavor structure of the meson. 
A'^ = A^^t"" denotes the gluon field, and the t° are the generators of the colour group. For 
large internal momenta the main contribution comes from those operators for which the 
quark fields and the path-ordered phase factor lie on a light-like straight line (xi — ^2)^ = 0. 
Such operators are known from the Operator Product Expansions [0]. They are highly 
singular and posses new properties. For simplicity let us write the corresponding Light-Cone 
Expansion for a toy model with scalar fields (p and scalar currents j: 

Tj{x)j{0) = J2 Cnun20n,,n2 + •••• , (3) 
ni,n2 

where the Cn^,„2(a;^) are singular coefficient functions, and Oni,n2 are the operators 

Or,,,n2 =■■ [{xd,r<P{xmxdX'<P{x)]\.=0 : , (4) 

and X is a light-like vector that approaches the vector x. This local light-cone expansion can 
be summed up to a nonlocal integral representation 

Tj{x)j{0)=f dKidK2F{x^,Ki,K2)0{Ki,K2) + (5) 

Jo 

where 



0(ki,K2) = : (j){Kix)(f){K2x) : (6) 
n7%2 ni\n2\ 



and 



1 

Cni,n2{x^) = i T / dKidK2F{x'^,Ki,K2)K'l^K2^. 

rii Ino! JO 



ni\n2 

Such a nonlocal light-ray operator product expansion has been independently introduced by 
and from different points of view. The operators 0(ki, 1^2) reflect the renormalization 
properties of all summed up local operators. The support of this operator is a light-like 
straight line which introduces new singularities and is responsible for the changed renormal- 
ization properties. The operator K2) depends on two parameters ki and K2 so that in a 



multiplicative renormalization scheme the Z-factors and the anomalous dimensions depend 
on two initial and two final parameters, i.e., 

0(ki, Ka)""" = dK\dK'^Z{Ki, K2] 4, 4)0(4, 4)""", (9) 







7(ki,/€2;4,K2) = (10) 

\ lo ('"^^'^^i' i^ii 4') j ^~^(4> 4'; 4> 4) 5 

where denotes the renormalization point. It is possible to derive the standard local anoma- 
lous dimensions from the nonlocal anomalous dimension. 

In a very formal way we may write the renormalization group equation for the operator 
0(ki, K2)'"^"(/i) = RTO{ki, K2)exp{iSj} according to 

I^^RTOMe'^' = (11) 



where R denotes the renormalization procedure, T the time ordering, and 5*/ is the interacting 
part of the action. The new variables are k± = (^2 ± ki)/2. For convenience we split the 
anomalous dimension of the operator into the anomalous dimension of the IPI vertex function 
(with two external momenta) 7(/i, k^) = 7(«+, ^ P^^^ which is proportional 

to the anomalous dimension of the quark field 7^. As short notations we use d'^isL = dn'j^dK.'_, 
and 5(2)(k - 4) = 6{k+ - k'^)5{k^ - k'_). 

2 Generalized Distribution Functions and Evolution 
Equations 

Before defining generalized distribution functions we repeat the definitions of the well-known 
meson wave function and the parton distribution function in QCD. 

The meson wave function ^""{x, Q^) which depends on the distribution parameter x (0 < 
X < 1) can be defined as the expectation value of a nonlocal (light-ray) operator lying on 
the light-cone § |, |]: 

$'^(x = i±^,g^)= (2.1) 

rdnJjiPl ^,,_(-p), ^ o|i?TO"(K_;n)e*^^|P > \,^=qK 
J 2'K[nP) 

Here \P > denotes the one-particle state of a scalar meson of momentum P, and 

0"(k_; n) = ■.^{-K_n){n-i)X'U{-K_h, K_n)i){K_n) : (2.2) 

is the light-ray operator with the same fiavor content. As the renormalization point we choose 
the typical large momentum Q of the basic process to which the wave function contributes. 
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The introduced factor l/(nP) compensates the n-dependence of the factor n7 in (2^). The 
vector n with = defines the hght-ray pointing in the direction of the large momentum 
flow of the process. The path ordered phase factor is taken along the straight line with 
direction n, and A" is a generator of the flavor group corresponding to the considered meson. 

Analogously we define the quark distribution function q°-{z,Q'^) [0 with the distribution 
parameter z. For simplicity we consider the flavor nonsinglet distribution only defined by 

q-{z, Q2) = / ^^-^ g2«_(nP). ^ p\RTO%ti^; n)e'^'\P > \^.^=Q^. (2.3) 
J 2T[{nF) 

If we choose the index a so that the matrix A'^ is diagonal, then for a positive resp. a 
negative distribution parameter z this function represents a linear combination of the quark 
resp. antiquark distribution functions 0. 



nP_ 




pO 


_pii 


nP+ 




PI 


-Pi 



The functions in ( p.lD and (|2.3|) are used for different physical processes and have dif- 
ferent interpretations. Nevertheless it seems to be natural to introduce the more general 
distribution amplitude 

q\t,r,l^')= (2.4) 
2n{nP^) 

with P± = P2± Pi. This function depends on the distribution parameter t and the quotient 
T = nP^/nP^ of the projection of momenta onto a light-like direction n, the renormalization 
point /i, and the scalar products of the external moments PiPj- For physical states \Pi > 
and \P2 > the additional variable r is restricted by 

« nP , , 

\n\ 

With the help of the a-representation it is possible to investigate the support properties 
of the function g"(t, r, /x^) with respect to the variable t. It turns out that q"'{t,T,Q'^) = 
for \t\ > 1. This distribution amplitude is of relevance for e.g., the description of 
nonforward processes near the forward case, where the parton picture with two distribution 
functions cannot be applied. On the other hand, it is possible to apply such a function in the 
limit P2 ^ too. Of course, then this state is not the vacuum state, but for the mathematical 
aspects of the connection of evolution equations for forward and nonforward processes this is 
very useful. In the limit of forward scattering we reach the "forward distribution amplitude" . 
This amplitude is real for Pi = P2, and for this reason it simultaneously plays the role of the 
parton distribution function in deep inelastic scattering (For the virtual Compton scattering 
amplitude approximated by a light-cone expansion the formation of its absorptive part is 
nontrivial for the hard scattering part only.). So we can perform two limits: 

$'^(x = (l + t)/2,Q2) = hm q%t,T,Q^) meson wave function, (2-6) 

T— » — 1 

q°'{z,Q'^) = \im q"'{z,T,Q'^) quark distribution function. 

Next we need evolution equations valid for the distribution amplitude ( |2.4| ). As input we 
can use the renormalization group equations ( |Tl| ) for the light-ray operators ( p. 21 ). As special 
cases we should find the Brodsky-Lepage (BL) and the Altarelli-Parisi (AP) equations. 
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For the derivation of the evolution equation we differentiate the general distribution 
amplitude q'^{t,T,fi^) with respect to the renormalization parameter. Thereby we take into 
account its representation in terms of matrix elements of light-ray operator ( |2.4| ). The 
differentiation of this operator can be performed with the help of its renormalization group 
equation (p!T|). A straightforward calculation starts from 

^±q'^(t,T,i^')= (2.7) 
djj 

r rf/t (nP+) ^^._inP^)t d ^ p^i^y^a^^^ ^ Q , K.' n)e''^ \P, > |nP_=.(nP+). 

J 2TT[nF^) dfi 
We obtain the evolution equation 

with the evolution kernel 

7(t,t') = Jdw--f{w+ = t'w--t,w-), (2.9) 

w+ = -t w„ = — , 2.10 

where 7(w+,w_) is directly related to the original anomalous dimension 7(/i, k^). Using 
translation and scale invariance we get 

7(/t+, /€_; k'-) = —7(0, 1, — , — ) = — 7(^+, ^-)- (2.11) 



This property is essential for reducing the two-variable renormalization group equation (|TI]) 
to a one- variable evolution equation ( |2.(j| ). From the renormalization group invariance of 
/ dt q"'{t, r, /i^) it follows that 



rft7(t,t') = 27^, (2.12) 
so that the standard -definition for the generalized function 

[7(t, t')]+ = lit, t') - 5(t - t') f dt" 7(t", t') = 7(t, t') - 27^(5(t - t') (2.13) 



arises in a very natural way. We see that instead of the original anomalous dimension there 
appears an evolution kernel that contains the original anomalous dimension as an essential 
input. It seems to be natural to denote the general evolution kernel as an extended BL-kernel. 
The reason for this is the following: restricting 7(t, t') to the parameter region \t\, \t'\ < 1, it 
coincides with the evolution kernel for the hadron wave function, but eqs. (|2.9| ) and ( p.ll|) 
provide us with a meaningful definition outside this region. 

In this way, we obtain evolution equations for all forward and nonforward matrix ele- 
ments. To our knowledge, in the literature up to now only two examples have been investi- 
gated, the case of forward scattering and the case of the meson wave function. Here, both 
cases are contained as limits P2 —>■ Pi or P2 —>■ 0. We write these limits down in a short 
form: 
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Evolution equation for the quark distribution function (forward scattering) with the 
AP-kernel P{z/z') 



(2.15) 



r^o 2r 



z z' 
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r r 



Evolution equation for meson wave functions with the BL-kernel VBL{x,y) 



$"(x = (l + t)/2,Q^) = lim g"(t, r,Q^), 

r— > — 1 

g'^$'^(^,Q') = /^'rfi/rB,.(x,i/;a.(Q2))$'^(l/,Q2). (2.16) 
VBL{x,y) = [7(2x-l,2y-l)]+|o<a=,y<i, (2.17) 

where QsiQ"^) denotes the running coupling constant. Here, the variables x = (1 + t)/2 
and y = {1 + t')/2 are restricted to < x,y < 1. Note, however, that in fact the 
quantum numbers must change discontinuously, so that the limit is only formal. 

Note that both kernels ( |2.15| ) and ( p.l7| ) are calculated separately in QCD. From the above 



considerations it is obvious, however, that they have a common origin, namely the anomalous 
dimension 7(w+, W-). 

3 The Extended BL-Kernel, Relations between the 
BL- and AP-Kernels 

We have obtained a new evolution kernel (extended BL-kernel) for distribution amplitudes 
containing the (restricted) BL-kernel and the AP-kernel as special cases. The standard 
perturbative calculations of the BL-kernel define it as an evolution kernel for a restricted 
range of the variables. It turns out that this kernel already contains the essential information 
and that an explicit continuation procedure can be prescribed. 

The first question is: What is the region in the (t,t')-plane where the kernel 'j{t,t') is 
defined? After some algebra we obtain the following representation in the (t,t')-plane [§]: 

7(t,0 = [e{t-tyii-t)-eit'-t)eit-i)]fit,t') (3.1) 

+ [9it' - t)9{l +t)- 9{t - t')e{-t - l)]/(-t, -t') 

+\e{-t - t')0{i + t)- e{t' + t)0{-t - i)]g{-t, t') 
+[e{t + t')e{i -t)- e{-t' - t)e{t - i)]g{t, -t'), 

where the functions f{t,t') and g{t,t') are given by 



i-t' 



f{t,t') = / dw- jIw^ = t — w_t' ,w_), 
Jo 

g{-t,t') = r ' dw_-f{w+ = -t-w_t',-w_). (3.2) 
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Note, that for \t\, \t'\ > 1, because of support restrictions, the lower boundary in the in- 
tegral representations (|3.2|) is not attained. However, in these regions only the differences 
f{t, t') — f{—t, —t') and g{—t, t') — g{t, —t') appear, so that such undetermined contribu- 
tions eliminate each other. It is clear, therefore, that the kernel is defined in the complete 
(t,t')-plane as it can be seen in Fig. 1. 

In one-loop approximation the extended BL-kernel reads 

V^l\x,y) = ^o{t = 2x-l,t' = 2y-l) = -^CF[Vo{x,yU (3.3) 

Voix,y) = e(l--]e(-]signiy)^(l-^]+{''^'!_]. 

\ yj \y/ y \ x-yj [y-^y j 

{cp is the known group theoretic constant). This result can be obtained by an extension 
procedure or by direct calculation of 7o(t,t'). 

As an interesting example we shall determine the AP-kernel from the known extended 
BL-kernel. According to eq. ( [^.15| ) it is just this new region {\t\, \t'\ > 1) which is essential 
for the determination of the AP-kernel. We first discuss the ^-structure and then consider 
the limiting process for the coefficient functions contained in eq. ( p.3|) . The limit r 0, 
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with X = z/t and y = 1/r, for the 5-f unctions are 




^ 9{z)9{l-z), (3.4) 

^ e{-z)e{i + z). (3.5) 

Related 6'-structures (structures that turn into each other under x = 1 — x and y = 
l~y) have the same hmit. Putting it all together we obtain an expression for the AP-kernel 
which follows from the extended BL-kernel ( |2.15D , ( |3.3| ), and 



and 



P{z) = -^c^[Po{z)U (3.6) 

ZTT 



Po{z) = limJr|-Vo(-,-) (3.7) 



+0' ' \T T 

1 + Z^ 



9{z)eii 



z 



The same procedure can be performed for the two-loop contributions to the BL-kernel. The 
calculations are much more complicated. The result obtained for the AP-kernel is equiva- 
lent to that of IP]. It also includes in a natural way the second order contributions stemming 
from the internal anti-quark lines. This constitutes a consistency check of the calculation of 
the diagonal local anomalous dimensions. The nondiagonal elements of the anomalous di- 



mension matrix have been controlled by exploiting conformal symmetry breaking |jTO|]. Note, 
that our regularization prescription for the extended BL-kernel, after performing the limiting 
process, turns automatically into the well accepted regularization prescription for the AP- 
kernel. 



4 Evolution Equations and Structure Functions for 
Two-Photon Processes 

Here we consider the virtual nonforward Compton scattering amplitude of a spinless particle 
(near forward direction) at a large off-shell behavior of the incoming photon. The considered 
region is a generalized Bjorken region, where (as in the case of deep inelastic scattering) this 
process is dominated by contributions from the light-cone. The kinematics of the process 
are the following: 

7*(gi) + H(pi) = 7*(g2) + H(p2), (4.1) 

P = Pi+P2 = {2E„6), q={l/2)iq^ + q2) (4.2) 

A = (p2-pi) = (gi-g2) = (0,-2p-), Ep = ^im^ + fP). 
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The last notations are the values of the momenta in the Breit frame. The generalized Bjorken 
region is given by 

u = Pq = 2Epqo ^ oo, = -q^ oo, (4.3) 

with the scaling variables 

^ = '^=-^ = ^ = —cos(j). 4.4 

Fq Fq Iv Ep 

To understand this process better we introduce the angle between the vectors p and q in the 
Breit frame by cos0 = pq/{\p\\(l\)- In terms of these variables we get 

= - - J| COS 0^ U, g2 = - + COS 0^ U. (4.5) 

In contrast to deep inelastic scattering the variable ^ is not restricted to < < 1. For 
example ^'2 = demands = —{\p\/Ep) cos (p. It may be shown that in this region the helicity 
amplitudes T(A', A) = e2{X')Tfj_i, x e'({X) are given by (l/2)£2(A')£:i(A)T^ for the transverse 
helicities and vanish otherwise. Therefore, only the trace of the scattering amplitude 

T,,{P,A,q)=zJd'xe'^^ < P^IT (j^ (0 J, |Pi > (4.6) 

has to be considered; {J^{x) = (1/2) : ip{x)'j^{\^ — \^ /\^)'ip{x) : is the electromagnetic 
current of the hadrons (for flavour SU{3)). 

In our special case the light-ray operator product expansion contains in leading order 
only the following quark operator 

0"'{x, Ki, K2) =: tp{Kix){x'y)U{KiX, K2x)'iIj{k2S:) : . (4.7) 

The singular coefficient functions Fa are determined perturbatively; with the Born approxi- 
mation of the coefficient functions we obtain finally 

T;(P, A,g) ^ 2 £dt - eaqUt,V;f^' = Q'), (4.8) 

where Ca = {2/9)Sao + (l/6)'^a3 + (l/6"\/3)'^a8; and the distribution amplitudes 

g^(t,r;/i2) = /^^e^"^'^* < P2|0'^(x, = 0),.\F, >, ^^^^ (4.9) 



contain the long range behaviour of the process (which is related to the hadron states |Pj >). 
This amplitude generalizes the parton distribution function of deep inelastic scattering. It 
has the following properties: 

• the support of g^(t, r, Q^) is restricted by \t\ < 1 and |r| < 1, 

• g"(t, Q"^) = qc(t, r = 0, Q'^) is essentially the parton distribution function used in deep 
inelastic scattering, 

• g^(t, r = 0, Q^) is real. 

This generalized distribution function satisfies our new evolution equation Eq. (|2.8|). This 



conference contribution collects results of the article 1 1 
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